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PREFACE 


This  technical  report  covers  the  work  performed  under  ONR 
Contract  N00014-80-C-0700  from  16  June  1982  to  15  June  1983  and 
is  the  third  report  published  under  the  program.  Dr.  Robert 
Whitehead  is  the  Scientific  Officer  for  the  Office  of  Naval 
Research.  Other  sponsors  are:  Naval  Air  Systems  Command,  the 
Army  Missile  Command,  the  Air  Force  Armament  Test  Laboratory, 
the  Air  Force  Flight  Dynamics  Laboratory,  NASA/Langley  Research 
Center,  and  NASA/Ames  Research  Center.  Mr.  Wallace  Sawyer  and 
Mr.  David  Shaw  of  Langley  are  the  liaison  officers  for  NASA. 

Dr.  Leon  Schindel  of  the  Naval  Surface  Weapons  Center  at  White 
Oak,  MD  is  the  Technical  Monitor. 
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1.  INTRODUCTION 


Under  Office  of  Naval  Research  Contract  No.  N00014-80-C- 
0700,  Nielsen  Engineering  &  Research,  Inc.  (NEAR)  has  been 
funded  to  obtain  the  necessary  fin  data  base  for  broadly  ap¬ 
plicable  engineering  prediction  programs  for  calculating  the 
aerodynamic  characteristics  of  body-tail  and  canard  (wing)- 
body-tail  missiles.  Under  the  same  contract,  the  data  base 
will  be  incorporated  into  one  specific  computer  program  called 
PROGRAM  MISSILE  (refs.  1  and  2)  which  is  valid  for  angles  of 
attack  up  to  45°  and  arbitrary  roll  angle.  The  first  year's 
work  involved:  (1)  selection  of  the  test  model  design,  test 
parameters  and  testing  sequence,  (2)  preliminary  investigation 
of  the  optimum  approach  for  data  handling,  i.e.,  preparing  the 
data  for  and  incorporating  it  into  PROGRAM  MISSILE  and  (3)  re¬ 
vising  the  equivalent  angle-of-attack  formulation  to  incorpo¬ 
rate  the  new  fin  deflection  data  base.  The  results  of  that 
work  are  described  in  the  first  year's  report  (ref.  3).  The 
second  year's  work  consisted  of  (1)  support  of  the  ongoing 
wind  tunnel  tests,  (2)  preparation  for  processing  the  data  to 
be  incorporated  into  the  data  base,  (3)  continued  improvement 
of  the  methods  used  in  MISSILE  and  (4)  continued  code  devel¬ 
opment.  The  results  of  that  work  are  described  in  the  second 
year's  report  (ref.  4).  The  third  year's  work  continued  the 
activities  of  the  second  year. 

Six  tunnel  entries  have  been  planned.  The  three  high 
Mach  number  (2.5  -  4.5)  tests  were  to  be  conducted  in  the  NASA/ 
Langley  Unitary  Plan  Tunnel,  Section  2.  The  first  and  second 
Langley  entries  were  completed  in  May  1982  and  July  1983  re¬ 
spectively.  The  three  low  and  intermediate  Mach  number  (0.6  - 
2.0)  tests  will  be  conducted  in  the  6-  by  6-foot  Supersonic 
Tunnel  at  NASA/Ames.  The  first  Ames  entry  was  completed  in 
January  1983.  The  second  Ames  entry  is  now  scheduled  for 
January  1984. 

At  the  writing  of  this  report,  all  of  the  6=0  data  for 
the  control  fins  have  been  obtained.  Because  of  problems  with 
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the  low-aspect-ratio  fin  balances,  no  data  for  those  fins  has 
been  obtained.  The  low-2R  fins  will  be  tested  in  the  third 
entry  at  NASA/Langley  and  in  the  second  entry  at  NASA/Ames. 
The  status  of  the  control  (6  ?  0)  data  is  given  in  the  table 
below.  An  "X"  means  that  the  data  have  been  obtained. 


Required  Control  Fin  Data  Base  (6  j-  0) 


\  M 

\  CJO 

FIN  N. 

Angle-of- 

attack 

range 

-  NASA/ARC  - ► 

< -  NASA/LRC  - »• 

o 

• 

00 

1.2 

2.0 

3.0 

n 

42 

S*=l 

X=l/2 

low 

X 

X 

high 

51 

PR=2 

X  =  0 

low 

X 

X 

X 

X 

X 

high 

X 

X 

52 

PR=  2 

X=l/2 

low 

X 

X 

X 

X 

X 

high 

X 

X 

53 

m=  2 

X=1 

low 

X 

X 

X 

high 

X 

X 

62 

m=  4 
X=2 


low 

high 


The  next  section  presents  some  representative  data  from 
the  first  entries  at  NASA/ARC  and  NASA/LRC.  In  the  third  sec¬ 
tion  an  analysis  is  presented  of  some  of  the  high  Mach  number 
control  effectiveness  data.  The  analysis  leads  to  some  impor¬ 
tant  conclusions  regarding  the  rational  modeling  to  be  used  in 
the  data-base  prediction  program  (MISSILE  3) .  In  section  4  a 
discussion  of  some  of  the  data  processing  methodology  is  pre¬ 
sented  including  changes  which  are  required  as  a  result  of  the 
findings  of  section  3.  In  section  5  an  extensive  analysis  of 
the  proposed  body  alone  methodology  is  presented.  Concluding 
remarks  are  given  in  section  6. 
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2.  REVIEW  OF  THE  DATA 


2.1  Control-Effectiveness  for  Transonic  Flow 

Some  sample  transonic  results  from  the  first  entry  at  NASA/ 
ARC  are  given  in  figure  1.  The  sign  conventions  are  given  in 
figure  2.  The  data  are  for  the  aspect  ratio  2  control  fins  at 
ac  =  20°  with  taper  ratios  equal  to  0,  1/2  and  1.  Only  the  nor¬ 
mal-force  coefficient  for  the  one  deflected  fin  is  shown.  Im¬ 
portant  things  to  note  are  (1)  the  change  in  taper  ratio  effects 
with  Mach  number  and  (2)  the  significantly  smaller  control  ef¬ 
fectiveness  for  5  =  40°  compared  to  6  =  -40°. 

2.2  Control-Effectiveness  for  Supersonic  Flow 

Typical  normal-force  coefficient  data  from  both  of  the 
entries  are  presented  as  a  function  of  bank  angle  in  figure  3 
for  the  deflected  fin  for  a  body  incidence  angle  of  20°.  The 
reference  area  is  the  fin  planform  area.  There  are  two  fea¬ 
tures  which  should  be  noted.  First,  the  shape  of  the  62  =  0 
curves  is  familiar  from  previous  work  (e.g.,  ref.  5).  The  max¬ 
imum  value  of  the  normal-force  coefficient  is  displaced  toward 
the  windward  side  from  the  horizontal  symmetry  position  of 
4»  =  0  to  roughly  0  =  2  5°  in  accordance  with  slender-body  theory, 
and  the  peak  values  of  decrease  with  M^  as  would  be  ex¬ 

pected.  Note  also  that  the  negative  values  of  CnF(B)  for  the 
<S2  =  0  cases  are  much  smaller  for  =  3.0  and  4.5  indicating  a 
change  in  vortical  flow  properties. 

The  second  feature  of  note  is  the  changing  control  effec¬ 
tiveness  with  increasing  Mach  number.  In  particular,  compare 
the  normal-force  coefficients  at  0  =  -90°  (leeward  plane) 
and  4>  =  +90°  (windward  plane).  For  S2  =  0,  CjgFjBj  at 
<J>  =  +90°  is  zero;  hence,  the  values  shown  for  S2  /  0  repre¬ 
sent  control  effectiveness  only.  For  M^  =  2.00,  the  control 
effectiveness  for  <J>  =  90°  is  slightly  greater  than  for  <J>  =  -90°. 

As  M  is  increased  the  control  effectiveness  at  <J>  =  -90°  decreases 
00 
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snarply  to  negligible  values  while  the  effectiveness  at  <p  =  90° 
increases.  Note  also  that,  for  certain  values  of  $  on  the  wind¬ 
ward  side,  is  greater  than  the  maximum  value  of  1.7  gi¬ 

ven  by  Hoerner  (ref.  6)  for  a  flat  plate  normal  to  the  free 
stream!  An  additional  set  of  data  for  the  same  fin  but  with 
ac  =  35°  and  Mb  =  3.0  is  given  in  figure  4.  Note  the  similarity 
of  the  curves  of  figure  4  with  those  of  figure  3(c)  which  have 
roughly  the  same  cross  flow  Mach  number. 

An  interesting  perspective  on  the  phenomenon  described 
above  can  be  gained  by  plotting  the  fin  normal-force  coefficient 
for  the  $  =  +90°  positions  as  a  function  of  cross  flow  Mach  num¬ 
ber,  M„  =  M  since  .  This  is  done  for  M  =  2.0,  3.0  and  4.5  and 
fin  deflection  angles  of  +20°  in  figure  5.  The  fin  normal-force 
coefficient  for  a  given  ac  is  normalized  by  its  value  at  otc  =  0 
to  give  the  curves  a  common  value  at  Mc  =  0.  Figure  5  shows 
that  the  data  correlate  fairly  well  with  cross  flow  Mach  number. 

It  is  obvious  from  figure  5  that  a  yaw  command  in  the  plus 
configuration  for  a  high  cross  flow  Mach  number  will  lead  to  an 
induced  rolling  moment  which  will  have  to  be  answered  by  the 
horizontal  fins.  However,  when  a  high  performance  missile  is 
pulling  a  high-g  maneuver,  much  of  the  available  horizontal  fin 
deflection  capability  will  be  used  to  trim  the  vehicle  in  the 
pitch  plane.  Hence,  little  will  be  left  to  counter  any  induced 
rolling  moment  due  to  a  yaw  command.  Clearly,  this  would  put  a 
limit  on  the  ac  which  the  vehicle  could  safely  reach  for  high 
cross  flow  Mach  numbers. 
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3.  ANALYSIS  OF  THE  CONTROL-EFFECTIVENESS  DATA 
FOR  SUPERSONIC  FLOW 


3.1  Correlation  of  the  Data 


The  large  variation  in  with  fin  bank  angle  and  the 

correlation  of  figure  5  prompts  us  to  consider  the  local  flow 
field  property  variation  around  the  body  alone .  To  do  this  we 
used  the  SWINT  marching  Euler  solver  (ref.  7)  which  has  the 
capability  of  shedding  and  convecting  vorticity  at  prescribed 
separation  lines  on  the  body.  The  separation  line  locations 
were  determined  from  a  correlation  by  Nielsen  (ref.  8)  made  from 
Landrum's  oil-flow  data  (ref.  9).  A  series  of  runs  were  made 
on  a  3-caliber  cone-cylinder  for  various  Mach  numbers  and  angles 
of  attack.  After  the  runs  were  made  a  postprocessor  was  used 
to  compute  the  local  dynamic  pressure,  q^,  and  Mach  number,  M^ . 

For  the  purposes  of  the  present  work,  we  will  use  values 
of  q^  and  M^  which  have  been  averaged  over  the  exposed  span  of 
the  fin,*  i.e. 


(1) 


^m 


m£  h 


s  -a 
m 


M^  ( 4> ,  r)  dr 


(2) 


Hence,  for  the  data  presented  here  the  integration  interval  was 
from  a  to  2a.  The  computed  results  for  q^  and  M^,  at  10  diameters 
aft  of  the  nose  tip  for  a  body  incidence  angle  of  20°  are  given 
in  figure  6. 

It  is  immediately  apparent  from  the  results  shown  in  figure 
6  that  the  variations  of  both  the  local  Mach  number  and  dynamic 


*The  averaging  was  applied  along  the  radial  line  defined  by  the 
exposed  fin  span  for  6  =  0  as  if  the  fin  were  immersed  in  the 
body  alone  flow  field. 


'.**  r-  t 
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pressure  contribute  to  the  behavior  displayed  in  figures  3-5. 
At  Ma  =  1.60,  for  example,  both  q^  and  M^  vary  only  slightly 
with  <J>  from  their  free  stream  values.  However,  as  M  increases, 
the  variation  with  <p  of  both  quantities  increases.  At  M^  =  4.63 
and  0  =  90°,  for  instance,  q^/q^  =  2.07  and  M^,  =  3.5.  Both  of 
these  flow  field  values  would  contribute  to  higher-than-expected 
values  of  CNf^B)  when  the  fin  is  deflected,  the  lower  Mach  num¬ 
ber  giving  a  higher  wing-alone,  normal-force  coefficient  slope. 

Given  the  above  results,  it  seems  reasonable  to  conjecture 
that  the  wide  variation  in  cNf(b)  seen  -*-n  figures  3  -  5  is  due 
to  local  flow  field  property  variations  caused  by  the  presence 
of  the  body  alone.  This  conjecture  can  be  checked  by  using  the 
equivalent  angle-of -attack  concept  of  Appendix  A  of  the  first 
year's  report  (ref.  3)  and  the  local  flow  field  property  values 
of  figure  6.  Since  only  fin  2  of  the  model  was  deflected  we 
can  write  equation  (A-10)  of  reference  3  as 


or 


eq2 


=  ot__  + 
eq2 


A2262 


A 

•22 


(3) 

(4) 


If  the  conjecture  is  correct,  it  should  be  possible  to  use 
equation  (4)  to  collapse  the  6  and  dependence  of  the  data 
in  figures  3  and  4. 

The  steps  which  were  used  to  check  the  conjecture  for  a 
given  body  angle  of  attack  and  free  stream  Mach  number  are  as 
follows : 

1.  For  particular  values  of  $ ,  M^  and  ac,  determine 

from  the  data  for  <5 2  =  0,  +20°,  +40°. 

2.  Multiply  the  results  of  step  1  by  qx/q £  for  the  same 

values  of  <f>,  M*  and  ac  to  obtain  Cjgp(g)  referenced  to 
the  (computed)  local  dynamic  pressure. 

3.  Interpolate  in  the  wing-alone  data  of  reference  10 

for  wings  with  the  same  planform  at  the  (computed) 
local  Mach  number  to  obtain  values  of  aeg 
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corresponding  to  the  results  of  step  2.  Ignore  those 
data  points  which  would  require  extrapolation  of  the 
wing-alone  curves  beyond  a  =  60°.  The  wing-alone 
curves  (ref.  10)  are  given  in  figure  7. 

4.  Apply  equation  (2)  to  the  results  of  step  3. 

The  results  for  the  procedure  above  applied  to  the  data 
of  figures  3  and  4  are  given  in  figures  8  and  9.  It  is  clear 
that  the  correlation  succeeds  in  taking  out  the  $  dependence 
of  the  control  effectiveness  very  well  except  for  the  bank  angle 
region  approximately  +35°  from  the  leeward  plane.  The  discrep¬ 
ancy  near  the  leeward  plane  increases  with  M^  from  hardly  notice¬ 
able  at  =  2.0  to  more  than  a  factor  of  2  at  MK  =  4.5.  Hence, 
it  seems  safe  to  conclude  that  the  Euler  solver  of  reference  7 
as  implemented  for  the  solution  of  figure  6  gives  incorrect  local 
flow  field  properties  near  the  leeward  plane  for  cross  flow  Mach 
numbers  greater  than  one.  However,  because  the  loads  under  those 
conditions  are  relatively  small  for  the  a/sm  ratio  tested  (0.5), 
the  error  would  not  be  particularly  important.  For  smaller 
values  of  a/sm,  the  error  may  be  more  important. 

3.2  Implications  for  the  Equivalent  Angle-of-attack 
Concept  and  MISSILE  3 

From  the  discussion  in  the  subsection  above,  it  is  clear 
that  the  aeg  concept  as  outlined  in  Appendix  A  of  reference  3  and 
in  reference  5  will  not  accurately  predict  control  effectiveness 
for  cross  flow  Mach  numbers  greater  than  0.5  because  it  is  im¬ 
plied  in  the  concept  that  q^  and  do  not  deviate  from  the  free- 
stream  conditions.  However,  with  a  reasonable  amount  of  effort, 
the  method  can  be  extended  so  that  it  can  still  be  used  in  the 
MISSILE  3  code  to  be  developed  under  the  present  contract.  To 
extend  the  aeq  concept  to  cross  flow  Mach  numbers  greater  than 
0.5  a  table  of  q^/q^  and  M^  as  functions  of  Mk,  ac,  $ ,  a/sm  and 
x/D  must  be  provided.  For  conditions  under  which  the  velocity 
parallel  to  the  body  axis  is  supersonic  everywhere,  such  a  table 
can  be  obtained  relatively  inexpensively  using  a  marching  Euler 
code  with  body  vortex  separation  (e.g.,  ref.  7). 
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The  region  in  which  a  marching  code  can  be  used  is  given 
in  figure  10.  As  part  of  the  work  for  the  present  effort,  an 
extensive  series  of  SWINT  (ref.  7)  runs  were  made  to  give  the 
necessary  body  vortex  properties  for  a  body  alone.  These  same 
runs  will  be  used  to  make  the  needed  MISSILE  3  table  for  q^/q^ 
and  M^.  For  the  region  bounded  by  Mc  =  0.5  and  the  shock- 
detachment  boundary,  an  unsteady  Navier-Stokes  code  will  be 
needed.  For  the  Mc  <  0.5  region  simple  point-vortex  "cloud" 
theory  will  be  used  to  obtain  the  body  vortex  properties  needed 
for  MISSILE  3  since  q^  and  M^  do  not  vary  significantly  from 
the  free  stream  values. 

The  actual  steps  required  to  extrapolate  from  the  fin-on- 
body  data  base  of  MISSILE  to  different  values  of  a/s  ,  x/D  and/ 
or  (Act)y^  are  as  follows: 

1.  For  the  M^,  ac ,  4>,  6  combination  of  interest,  deter¬ 

mine  Cn_,dv  from  data  base. 
t  (B) 

2.  Determine  q^/q^  and  M^  for  the  <J>  location  of  interest 

for  the  a/sn  and  x/D  of  the  fin  in  the  data  base. 

3.  Divide  the  CNp(B\  of  step  1  by  the  q^/q^  from  step  2 

to  normalize  it  properly  for  use  in  the  aeq  method. 

4.  For  the  wing-alone  curve  corresponding  to  the  M^  of 

step  2,  obtain  aeq  for  the  normal-force  coefficient 

resulting  from  step  3. 

5.  Using  equations  A{9)  and  A(10)  of  reference  3,  deter¬ 

mine  the  new  aeq  for  the  conditions  of  interest. 

6.  Determine  q^/q®  and  M^  for  the  <t  location  of  interest 

for  the  a/sm  and  x/D  of  interest. 

7.  For  the  aeq  of  step  5  determine  the  C^w  corresponding  to 

the  M^  of  step  6. 

8.  Multiply  the  of  step  7  by  the  q^/q^  of  step  6  to  get 

the  CnF(b)  of  interest  referenced  to  free  stream  con¬ 
ditions  . 

The  steps  above  are  illustrated  in  figure  11. 
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4. 


DATA  PROCESSING 


The  portion  of  the  data  processing  program  which  handles 
6=0  cases  without  vortex  asymmetry  has  been  checked  out  using 
the  data  for  fin  42  from  the  first  entry  at  NASA/LRC •  An  example 
of  the  tabular  output  for  fin  normal-force  coefficient  is  given 
in  table  1.  Results  from  the  graphics  routine  which  is  used  to 
check  the  normal-force  coefficient  tables  (see  subsection  3.4  of 
reference  4)  are  given  in  figure  12.  Note  the  more  "peaky"  dis¬ 
tribution  for  the  higher  Mach  numbers. 

Another  task  of  the  data  processing  program  is  the  corre¬ 
lation  of  fin  center-of-pressure  location  with  fin  normal-force 
coefficient.  The  algorithm  used  is  a  piecewise  linear  curve  fit 
to  the  windward  fin  data  (to  minimize  vortex  effects)  with  smoothing 
based  on  minimization  of  the  second  derivatives.  Typical  results 
are  given  in  figure  13.  Note  that  the  standard  deviation  is  also 
given.  It  can  be  used  in  the  new  version  of  MISSILE  to  estimate 
the  hinge-  and  bending-moment  errors  due  to  spread  in  the  cor¬ 
relation.  It  was  found  that  the  amount  of  smoothing  did  not 
significantly  change  the  standard  deviations  of  the  correlations. 
The  quantity  "GAMMA"  is  the  relative  amount  of  smoothing  used. 

For  the  correlations  of  figure  13,  it  was  assumed  that  the 
measured  normal-force  coefficients  are  the  appropriate  values  to 
use  and  that  the  Mach  number  seen  by  the  fins  is  essentially  M^. 

In  the  previous  section  it  is  demonstrated  that  the  above  as¬ 
sumptions  do  not  hold  for  high  cross -flow  Mach  numbers. 
Consequently,  it  is  expected  that  the  correlation  procedure  used 
above  will  have  to  be  altered  so  that  nonlinear  flow  field  ef¬ 
fects  can  be  taken  into  account.  Preliminary  work  indicates  that 
the  correlations  will  change  significantly  and  that  the  standard 
deviations  of  the  correlations  will  be  reduced. 
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5.  ANALYSIS  OF  PROPOSED  BODY  ALONE  NORMAL-FORCE  AND 
PITCHING-MOMENT  METHODOLOGY 


In  the  following  we  extensively  check  two  different  formu¬ 
lations  of  the  simplest  crossflow  drag  method  for  predicting  the 
normal  force  and  location  of  the  center-of-pressure  for  commonly 
encountered  missile  bodies  alone  using  empirical  input  for  CNq  and 
the  center-of-pressure  location  for  small  a.  The  ranges  of  the 
parameters  for  pointed  tangent-ogive  noses  in  combination  with 
right-circular-cylindrical  bodies  for  which  the  formulations 
were  checked  are  as  follows: 

+ 


/ 


0.5  <  £n/D  <4  °.6  <  <  5 

0  <  lA/ D  <20  0°  <  a  <  45° 

5.1  Candidate  Formulations  -  Normal  Force  Coefficient 
The  common  starting  point  for  the  formulations  considered 

is* 


CN  = 


N 


q-A 


=  sin2cx  +  nc 


n  r- 

oo  ref 


d  A  , 
c  ref 


.  2 
sin  a 


(5) 


*The  possible  advantage  of  using  an  axial  variation  of  the  cross 
flow  drag  coefficient  has  not  been  considered  in  the  present 
study. 
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where  N  is  the  resultant  normal  force  on  the  body,  cdc  is 
the  steady-state  crossflow  drag  coefficient  for  an  infinite 
circular  cylinder,  n  is  a  correction  factor  for  finite  cylinder 
length,  and  Ap  is  the  planform  area  of  the  body.  The  ref¬ 
erence  area  used  is  the  cross  sectional  area  of  the  cylindrical 
portion  of  the  body.  This  relation  is  a  straightforward  com¬ 
bination  of  slender-body  theory  (sin2a)  and  a  term  representing 
the  integrated  force  on  the  configuration  due  to  viscous  cross- 
flow  separation.  In  practice,  the  slender-body  term  has  often 
been  modified  by  the  factor  cos (a/2),  following  Ward  (ref.  11). 
With  this  modification,  the  above  relation  is  that  used  in  Ref¬ 
erence  12.  Note  that  at  a  =  90°, 

CN  =  ncdc  Ap/Aref  (6) 

In  the  present  work,  one  of  our  candidate  formulations 
(Method  A)  is  obtained  from  equation  (5)  as 


°Na  a  Ap 

CM  =  -  sin2a  cos  -  +  cH  — — 

N  2  2  dc  Aref 


.  2 
sm  a 


Inclusion  of  the  additional  factor  CN  /2  in  the  slender-body 
term  allows  for  empirical  input  of  the  slope  of  the  normal-force 
curve  at  a  =  0°,  instead  of  relying  on  the  slender-body  value 
(CNa  =  2).  Additionally,  we  have  set  n  =  1  to  be  used  in  con¬ 
junction  with  the  modified  curve  for  cdc  discussed  below. 

The  second  candidate  formulation  is  derived  from  Equation 
(5)  as  follows.  The  slender-body  contribution  to  CN  is  impor¬ 
tant  only  at  small  a,  for  which 

CN 

— — a  sin2a  cos  -  =  Cv,  sina  (8) 

2  2 

To  maintain  the  result  of  equation  (5)  while  using  the  approxi¬ 
mation  of  equation  (8) ,  we  must  write 
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CM  =  C.,  sina  +  c,  — *- - CM 

N  N  d  *  ,  N 

a  c  Aref  a 


2 

-  C„  sm  a 


Equation  (9)  will  be  called  Method  B. 


5.2  Candidate  Formulations  -  Center  of  Pressure 

Both  of  the  methods  to  be  presented  for  calculating  the 
location  of  the  center  of  pressure  (xCp)  are  based  on  the  as¬ 
sumption  that  the  contribution  to  CN  from  the  CN^  term  (Cnsb* 
the  term  linear  in  sin  a)  acts  for  all  a  at  the  center-of- 
pressure  location  of  the  configuration  at  zero  angle  of  at¬ 
tack  (x„_  ) ,  and  the  contribution  from  the  cross-flow  drag 
(CNv#  the  term  involving  sin  a)  acts  at  the  centroid  of  the 
planform  area  (xc) .  Thus, 


x 

C  +  —  C 

nsb  d  nv 


For  convenience,  the  alternative  formulations  are  sum¬ 
marized  in  the  table  below. 


CANDIDATE  METHODS 


C  +  C 

nsb  nv 


cp  X 

*0  C 

- —  C  +  —  c 

D  Nsb  D  LNv 


Method 


SB 


w 


V 


A 


B 


sin2a 


cos 


a 

2 


2 


sm  a 


C.,  sina 
N 

a 


.  2 
sm  a 


For  the  comparisons  with  data  the  methods  will  use  a  common 
data  base  for  CN^  and  xCpQ  (described  next)  and  the  same  rela¬ 
tion  for  c^^  (Figure  14) .  Figure  14  is  adapted  from  Reference  13 
with  the  irregular  behavior  of  c^  in  the  region  of  M  =  1  re¬ 
placed  with  the  faired  curve  shown.  For  the  comparisons  to  be 

shown,  the  branch  of  the  curve  for  turbulent  flow  for  M  <  .5 

c 

has  been  used. 


5.3  Data  Base  for  C.,  and  x„„ 

Not  CP0 

A  literature  survey  was  made  to  identify  existing  data  for 
these  quantities.  The  search  was  limited  to  data  for  bodies  with 
pointed  tangent-ogive  noses.  The  table  below  summarizes  the  results 
for  CNa-  The  situation  for  xCpo  is  very  similar.  Examination  of 
this  table  reveals  that  the  data  base  of  Reference  14  very  nearly 
covers  the  ranges  of  interest  of  £N/D,  £A/D  an<^  Ma>'  an^  that 
none  of  the  other  data  bases  extend  the  coverage  in  any  signifi¬ 
cant  way.  The  data  base  of  Reference  14  was  therefore  selected 
for  use  in  this  study.  To  assess  the  level  of  agreement  among 
the  data  bases  where  there  is  overlap,  a  series  of  comparisons 
were  made.  Figures  15  and  16  are  samples  of  those  comparisons 
for  CNa  and  xCpQ,  respectively,  for  the  data  of  References  13 
and  14.  On  the  basis  of  these  and  the  other  comparisons  not 
shown,  it  appears  that  the  data  from  the  different  sources  are 
generally  within  a  band  of  about  10%  for  CN  and  0.25  for  xCpQ/D. 


These  discrepancies  obviously  limit  the  ultimate  accuracy 
achievable  by  any  of  the  methods  being  considered. 


DATA  BASES  FOR  CN  FOR  BODIES  WITH  POINTED 

1Na 

TANGENT-OGIVE  NOSES 
Ranges  of  Parameters  Covered 


References 

VD 

VD 

M 

00 

13 

1.5,2. 5, 3. 5 

6-18* 

0.8-1. 2 

2. 5, 3. 0,3. 5, 4.0 

4-10 

1. 5-3.0 

14 

0.5, 1.5, 2. 5, 3. 5 

6-18 

0. 8-4.0 

15 

3,5,7 

0-10 

3-6.28 

16** 

2,3,4 

0 

2-4.5 

17 

2.84 

0-11 

1.36-4 

18 

1,2.4 

4,8 

0.6-1. 5 

5.4  Comparisons  with  Data 

To  evaluate  which  of  the  two  methods  previously  outlined 
is  most  accurate,  predictions  were  made  for  various  combinations 
of  parameters  for  bodies  with  pointed  tangent-ogive  noses  for  which 
there  are  experimental  data .  These  predictions  are  compared  to  the 
data  in  Figures  17-29.  Each  of  these  figures  is  for  a  particu¬ 
lar  combination  of  £N/D ,  £A/D  and  M^r  part  (a)  of  each  presents 
the  comparison  for  CN ,  part  (b)  for  xCp/D,  The  following  con¬ 
figurations  are  investigated: 


*For  £N/D  =  1.5,  data  exist  for  6  <_  £A/D  <  20. 

**Calculated  using  3-D  Method  of  Characteristics. 


Figure  No. 


VD 

VD 

1.5 

13 

3.0  3.667 

1  1 

3.5  7 


0.8 

17 

1.2 

18 

2.0 

19 

1.2 

20 

2.0 

21 

00 

• 

o 

22 

1.2 

23 

2.0 

24 

1.6 

25 

4.63 

26 

0.9 

27 

1.2 

28 

2.0 

29 

Examination  of  the  comparisons  for  CN  in  these  figures 
reveals  that  the  data  are  best  predicted  by  Method  B,  and  that 
the  discrepancy  between  that  method  and  the  data  is  generally 
less  than  10%.  The  only  exception  to  this  observation  is  for 
the  high-a  end  of  the  data  for  M^  =  0.8  in  Figure  17(a).  For 
this  case,  there  is  the  possibility  that  the  data  contain  the 
effects  of  asymmetric  vortex  shedding,  an  effect  clearly  not 
accounted  for  in  the  simple  methods  considered  here.  Note  that 
even  in  this  case,  Method  B  is  closer  to  the  data  than  Method  A. 

With  respect  to  the  location  of  the  center  of  pressure, 
the  agreement  with  data  of  the  predictions  is  considerably  more 
variable.  There  is  no  clear  advantage  for  any  of  the  methods 
over  the  others,  so  for  consistency,  Method  B  is  selected  for 
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1 

use  here  as  well  as  for  CN.  Figures  17  through  29  show  that 
the  prediction  of  Method  B  for  xCp  is  within  0.6D  of  the  meas¬ 
ured  location  for  nearly  all  of  the  cases  considered,  but  the 
discrepancy  is  up  to  two  to  three  times  that  level  in  the  cases 
of  Figures  17  and  18.  As  for  CN,  however,  in  these  cases  of 
lowest  accuracy,  Method  B  is  still  superior  to  Method  A. 

5.5  Conclusions 

The  following  formulations  (used  with  the  data  base  of 
Reference  14  and  the  cd  curve  of  Figure  14)  yield  predictions 
of  CN  and  xcp/D  that  are  generally  within  10%  and  0.6,  respec¬ 
tively,  of  measurements  over  the  ranges  of  parameters  shown: 


C 


N 


sina 


+ 


2 


sm  a 


(11) 


(CN  sina ) 
a 


+ 


.  2 
sm  a 


(12) 


0.5  <  £n/D  I  3-5 

6  £  ^a/d  -  18 

0 . 8  <  M  <4.0 
—  00  — 

0  <  a  <  45° 

These  ranges  are  set  by  the  data  base  of  Reference  14.  It  is 
estimated  that  the  small  extrapolations  necessary  to  cover  the 
complete  parameter  space  of  interest  can  probably  be  made  with¬ 
out  appreciable  loss  of  accuracy. 
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6 .  CONCLUDING  REMARKS 


At  the  end  of  the  third-year's  work,  the  following  remarks 
can  be  made: 

1.  For  M<*>  =  1.2,  CnF(B)  appears  to  be  a  fairly  strong 
function  of  A  for  0  <  A  <_  1/2,  but  it  is  only  a  weak  function 
of  A  for  1/2  £  A  £  1. 

2.  For  cross  flow  Mach  numbers  greater  than  0.5,  the  local 
properties  of  the  flow  field  about  a  body  have  a  very  strong  in¬ 
fluence  on  C^p(B)  resulting  in  the  possibility  of  large  control 
cross-coupling  effects. 

3.  The  variation  of  q^  and  M^  around  a  body  alone  for 
Mc  >  .5  will  have  to  be  taken  into  account  for  correlation  of 
center-of-pressure  data  for  the  data  base  and  for  extrapolating 
from  the  data  base  to  other  configurations  using  the  aeg  concept. 

4.  The  simplest  version  of  the  cross  flow  drag  method 
appears  to  be  capable  of  predicting  body-alone  normal  force  and 
center-of-pressure  location  to  within  10%  and  0.6  diameters 
respectively  provided  the  vortex  shedding  is  symmetric. 
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f  control  effectiveness  for  leeward 
dward  fin  positions. 


(a i  Averaue  local  dynamic  pressure. 
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aft  o 


fe.-  Computed  flow  field  properites  10  diameters 
f  the  nose  tip  of  a  3-cal iber-cone  cylinder 
at  20°  anqlo  of  attack;  a  /s^  =  0.5. 
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Figure  8.-  Concluded. 
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Figure  17.-  Comparison  of  predictions  and  data 
; n/D  =1.5,  -a/D  =  13,  =  0.8 
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Figure  27.-  Comparison  of  predictions  ar.d  dat 
■  iNj/D  =  3.3,  'a/D  =  7,  M,  =  0.0 
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NOMENCLATURE 


1 


PR 

Aref 

cdc 

cr 

CN 

CNF(B)  i 

CNSB 

CNv 

CNW 

CNU 

D 

fA 

Mc 

M^ 

M^ 


N 


V 

0O 


Body  Radius 
Planform  area 

Aspect  ratio  of  wing  alone 
Reference  area 
Cross  flow  drag  coefficient 
Fin  root  chord 
Normal-force  coefficient 

Normal-force  coefficient  of  fin  i  in  the  presence 
of  a  circular  body 

Linear  portion  of  the  normal  force  acting  on  a 
body 

Nonlinear  portion  of  the  normal  force  acting  on 
a  body 

Normal-force  coefficient  of  wing  alone 

Derivative  of  the  normal-force  coefficient  with 
respect  to  a 

Diameter  of  the  cylindrical  portion  of  a  body 
Length  of  the  cylindrical  portion  of  a  body 
Nose  length 

Cross  flow  Mach  number 

Mach  number  at  a  point  in  a  body-alone  flow  field 

Local  Mach  number  in  presence  of  body  alone  aver¬ 
aged  over  exposed  span  of  fin  (equation  (2)) 

Free-stream  Mach  number 

Normal  force  acting  on  the  body 

Dynamic  pressure  at  a  point  in  a  body-alone  flow 
field 

Local  dynamic  pressure  in  presence  of  body  alone 
averaged  over  exposed  span  of  fin  (equation  (1)) 

Dynamic  pressure  of  the  free  stream 

Radial  distance  from  body  axis 

Radius  of  curvature  of  tangent-ogive  nose 

Exposed  semispan  of  a  fin 

Semispan  of  fin 

Free-stream  velocity 
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NOMENCLATURE  (CONTINUED) 


X 

Distance  aft  of  nosetip 

xc 

Axial  location  of  the  centroid  of  the  planform 
area  of  a  body 

xcp 

Axial  location  of  the  center  of  pressure 

XCPo 

Axial  location  of  the  center  of  pressure  for  the 
linear  portion  of  the  normal  force 

^cp 

Lateral  location  of  the  center  of  pressure  of  a 
fin 

a 

Angle  of  attack 

ac 

Angle  between  body  axis  and  wind  velocity  vector 

aeqi 

Equivalent  angle  of  attack  of  fin  i;  i.e.,  angle 
of  attack  of  wing  alone  which  gives  same  normal- 
force  coefficient  as  that  of  fin  i 

<*eqi 

Equivalent  angle  of  attack  of  fin  i  if  all  fins 

are  undeflected 

*i 

Deflection  of  fin  i,  positive  when  the  leading 
edge  is  rotated  toward  the  leeward  side  of  the 
body 

( Aa) 

Average  angle  of  attack  induced  on  fin  i  by 
vortices 

Aaeq 

Increment  in  aeg 

n 

Length  factor  in  cross  flow  drag  theory 

X 

Fin  taper  ratio 

Aji 

Fin  deflection  factor 

♦i 

Roll  angle  of  fin  i 

